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Introduction 


Breast  cancer  spread  into  the  lymph  nodes  (LN)  is  the  single  most  deadly 
prognostic  factor  identified  to  date.  Understanding  factors  that  regulate  cancer¬ 
cell  colonization  of  the  lymphatics  is  essential  to  prevent  malignant  invasion. 

Most  breast  cancers  are  estrogen  (ER)  and/or  progesterone  receptor  (PR) 
positive,  and  the  majority  of  breast  cancers  that  spread  to  the  LN  maintain  ER 
and  PR  expression.  However,  due  to  lack  of  experimental  models  of  ER+  breast 
cancer  metastasis,  little  is  known  about  the  roles  of  hormones  in  breast  cancer 
spread  to,  and  growth  in,  LNs. 

Hypotheses:  I  postulate  that  breast  cancer  metastasis  to  LNs  is  dependent  on 
tumor  growth  and  tumor  size.  Second,  I  predict  that  tumor  spread  to  LNs  is 
dependent  on  tumor  associated  lymphangiogenesis.  Finally,  I  hypothesize  that 
ER  expression  and  transcriptional  function  is  maintained  in  tumors  and  LN 
metastases. 

Project  goals:  I  plan  to  develop  and  then  define  models  of  ER+  breast  cancer 
metastasis,  in  preparation  for  studies  of  hormone  action  that  will  address 
hormonal  regulation  of  key  genes  involved  in  lymphatic  development,  LN  spread, 
and  growth  of  cancer  cells  within  the  LNs. 

Approved  STATEMENT  OF  WORK: 

Task  1.  Develop  and  characterize  all  cell  lines  used  in  the  study  (Months  1- 

6). 

A.  Make  ZsGreen  and  DsRed-Express  polyclonal  stable  cell  lines. 

B.  Quantitation  of  ER  and  PR  protein  expression  in  each  lines. 

Task  2.  Assess  the  impact  of  systemic  estrogen  and  progesterone  on  ER+/- 
PR+/-  breast  cancer  lymphatic  spread  (Months  6-36). 

A.  Grow  and  determine  the  rate  of  lymphatic  invasion  using 
homogeneous  or  dual  colored  mixed  ER+/-,  PR+/-  tumors  with  no 
systemic  hormones,  E,  P,  E+P. 

B.  Perform  immunohistochemistry  on  tumors,  lymphatics,  and  LNs 
to  understand  effects  of  E  and  P  on  cancer  cells,  lymphatics,  and 
LNs. 

Task  3.  Determine  hormonal  influence  to  tumor  lymphangiogenesis 
(Months  8-36). 

A.  Determine  hormone  influence  on  gene  expression  of 
lymphangiogenic  factors. 

B.  Immunohistochemically  quantitate  lymph  vessels  in  each  tumor 
line  and  determine  the  influence  of  hormones  on  lymphatic  vessel 
development. 
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Body 


Task  1.  The  majority  of  the  work  completed  during  the  second  year  of  this  award 
is  summarized  in  my  2'^'^  and  3'^'^  first-author  publications  enclosed  within  the 
appendix  (1-2).  In  these  papers  I  describe  the  studies  that  complete  and  expand 
on  task  1  of  the  approved  statement  of  work.  This  includes;  (a)  using 
fluorescently  labeled  ER+  breast  cancer  cells  (MCF-7  and  T47D)  to 
immunohistochemically  quantify  the  amounts  of  ER  and  PR  expression  in  the 
breast  cancer  cells  grown  in  vivo  (see  appendix,  Harrell  et  al.,  Cancer  Research 
2007,  figure  2).  This  work  expands  on  my  last  year’s  first  author  publication  (3). 
By  using  estradiol  withdrawal  post  LN  metastasis  I  was  able  to  define  the  genes 
that  are  estrogen  regulated  in  primary  breast  tumors  compared  to  their  LN 
metastases  (see  appendix,  Harrell  et  al..  Cancer  Research  2007  figures  3-6).  We 
confirmed  that  ER  were  highly  expressed  in  mammary  gland  tumors  of  all 
estradiol-free  mice,  and  they  behaved  normally,  with  ligand  dependent 
downregulation  upon  estradiol  treatment  coupled  to  PR  upregulation.  However, 
these  expression  patterns  were  altered  in  metastases  within  the  lymphatic 
vessels  or  in  LNs,  with  poor  ER  downregulation  and  decreased  PR  induction 
(see  appendix,  Harrell  et  al..  Cancer  Research  2007  figure  2).  Thus,  the 
lymphatic  environment  makes  ER+  cancer  cells  estrogen  insensitive.  To  study 
this  further,  I  gene  expression  profiled  estrogen-treated  and  estrogen-withdrawn 
tumors  and  matched  LN  metastases.  For  my  first  study  (1)  the  microarray 
analysis  was  performed  on  laser-captured  MCF-7  cells  from  tumors  and  LN 
metastases,  instead  of  on  whole-tissue.  Interestingly,  the  array  data  confirmed 
the  immunohistochemistry  in  that  fewer  genes  were  E  regulated  in  LNs  as 
compared  to  primary  tumors  (see  appendix,  Harrell  et  al..  Cancer  Research  2007 
figures  3,  4,  6).  Furthermore,  some  genes  that  were  regulated  in  one  direction  by 
estradiol  in  the  primary  tumor  were  regulated  in  the  opposite  direction  by 
estradiol  in  the  LN  metastasis  (see  appendix,  Harrell  et  al..  Cancer  Research 
2007  figures  4,  6).  These  findings  are  very  important  as  they  provide  further 
evidence  that  the  LN  microenvironment  promotes  estrogen  resistance.  My 
second  study  from  this  2nd  year  of  funding  contrasted  how  whole-tissue  versus 
laser-captured  material  alters  LN  metastasis  gene  signatures  (2).  Interestingly,  I 
found  that  mouse  RNA,  which  contaminates  whole-organ  microarray  analyses, 
impacts  the  final  data  (see  appendix,  Harrell  et  al.,  ClinExpMet  2008  figures  2,  3). 
The  mouse  hybridizes  to  and/or  dilutes  out  the  human  RNA,  on  the  human 
specific  chips.  I  found  that  whole-tissue  and  laser-captured  isolated  RNA  yield 
totally  distinct  “metastasis  signatures”  and  I  speculate  that  the  laser-capture- 
derived  data  are  more  accurate.  These  studies  complete  task  1. 

Task  2.  The  majority  of  the  work  to  complete  task  2  is  currently  underway.  Over 
the  next  12  months  I  will  continue  studies  that  investigate  how  progesterone  (or 
medroxy-progesterone  acetate  (MPA))  alone  or  in  combination  with  estradiol 
affect  ER+  tumor  growth  and  metastasis.  The  findings  of  reduced  PR  expression 
(from  task  1  above)  in  metastases  suggests  that  treating  mice  with  estradiol  + 
progesterone  may  not  affect  cancer  cells  any  differently  once  they  become 
metastatic — as  the  amount  of  PR  that  can  bind  progesterone  is  diminished. 
Indeed,  just  a  few  weeks  into  these  studies,  I  am  not  observing  any  effect  of 
progesterone/MPA  on  tumor  growth  and/or  metastasis. 


5 


Task  3.  The  third  task  in  my  approved  statement  of  work  is  to  determine 
hormonal  influence  on  tumor  lymphangiogenesis.  In  my  year-one  first-author 
publication  I  show  that  estradiol  dependent  tumor  growth  correlates  with 
intratumoral  lymphangiogenesis  (reference  3,  figure  4);  intratumoral  lymphatic 
vessels  were  only  identified  in  estradiol  treated  MCF-7  tumors.  Peritumoral  LVs 
were  found  surrounding  nearly  all  tumors,  regardless  of  hormone  treatment,  and 
likely  derived  from  the  normal  lymphatic  vasculature  present  in  the  mammary 
gland.  To  further  investigate  hormonal  roles  in  lymphatic  vessel  development  I 
set  up  collaborations  with  Dr.  Kenneth  Korach  at  the  National  Institutes  of 
Environmental  Health  Sciences  and  Dr.  Darryl  Russell  at  the  University  of 
Adelaide  (South  Australia).  Through  RT-PCR  analyses  we  found  that  ER  alpha  is 
required  in  the  ovary  for  proper  lymphatic  remodeling  that  occurs  after  ovulation, 
whereas  ER  beta  has  no  role  in  this  process  (Figure  1).  This  is  the  first 
demonstration  that  estrogen  controls  lymphangiogenesis.  This  project  is  being 
expanded  upon  by  my  collaborators  and  we  expect  a  publication  within  the  next 
year.  These  studies  complete  task  3. 


Expression  of  the  lymphatic  marker  LYVE  after  ovulation 
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Genotype  and  Treatment 


Figure  1:  Postovulatory  lymphatic  vessel  remodeling  is  deficient  when 
ER  alpha  is  deficient.  Immature  WT,  ER  alpha  knockout,  and  ER  beta 
knockout  mice  were  stimulated  to  ovulate  and  ovaries  were  removed, 
prepared  into  RNA,  and  RT-PCR  was  performed  for  the  lymphatic 
marker  LYVE.  n=4  for  each  group  and  standard  deviation  is  shown. 


Ongoing  studies.  As  listed  above,  current  studies  are  underway  to  complete 
task  2. 
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Key  Research  Accomplishments 


•  Developed  a  unique  model  to  assess  estrogen-regulated  genes  in  tumors 
and  LN  metastases. 

•  Identified  that  the  LN  environment  alters  ER  expression  and  function. 

•  Performed  laser-capture  and  microarray  analyses  of  cancer  cells  from 
tumors  and  lymph  node  metastases  and  found  that  LN  metastases  are 
relatively  insensitive  to  estrogen  compared  to  primary  tumors. 

•  Identified  genes  that  are  regulated  by  the  microenvironment,  by 
hormones,  or  by  both. 

•  Identified  genes  that  are  regulated  by  estrogen  in  opposite  directions 
when  cancer  cells  are  in  the  mammary  gland  as  compared  to  LN 
metastases. 

•  Found  that  laser-captured,  as  compared  to  whole  tissue  extracted  RNA, 
yield  microarray  analyses  that  are  more  accurate.  These  findings  are  likely 
due  to  whole  tissue  samples  being  composed  of  different  amounts  of 
human  cancer  cells. 


Reportable  Outcomes 


1 - First  Author  Manuscripts 

Harrell  JC,  Dye  WW,  Harvell  DM,  Pinto  M,  Jedlicka  P,  Sartorius  CA,  Horwitz 
KB.  Estrogen  insensitivity  in  a  model  of  estrogen  receptor  positive  breast 
cancer  lymph  node  metastasis.  Cancer  Res.  2007  Nov  1;67(21):1 0582-91. 

Harrell  JC,  Dye  WW,  Harvell  DM,  Sartorius  CA,  Horwitz  KB.  Contaminating 
cells  alter  gene  signatures  in  whole  organ  versus  laser  capture 
microdissected  tumors:  a  comparison  of  experimental  breast  cancers  and 
their  lymph  node  metastases.  Clin  Exp  Metastasis.  2008;25(1)  81-8. 

2- Other  Manuscripts 

Christenson  KL,  Micalizzi  DS,  Coletta  RD,  Jedlicka  P,  Harrell  JC,  Horwitz  KB, 
Billheimer  D,  Heichman  K,  Welm  A,  Ford  HL.  The  developmental  regulator  Sixl 
promotes  an  epithelial  to  mesenchymal  transition,  induces  metastasis,  and  predicts 
poor  prognosis  in  many  human  cancers.  Submitted,  Nature  Medicine,  February 
2008. 

Horwitz  KB,  Dye  WW,  Harrell  JC,  Jedlicka  P,  Sartorius  CA.  Steroid  Receptor 
Negative  Breast  Cancer  Stem  Cells  are  Expanded  by  Progesterone  and  Generate 
Receptor  Positive  Tumor  Cell  Populations.  PNAS,  In  Press,  February  2008. 

3- Degree  Compietion  &  Future  Empioyment 

My  thesis  was  completed  in  the  summer  of  2007.  Upon  graduation  I 
stayed  in  the  Horwitz  lab  to  complete  this  project  and  to  set-up  my  postdoctoral 
lab.  I  have  since  accepted  an  offer  to  postdoc  with  Dr.  Charles  M.  Perou  at  the 
University  of  North  Carolina,  and  will  start  in  his  lab  towards  the  end  of  2008. 
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Conclusions 


The  second  year  of  the  funding  (March  1  2007-March  1,  2008)  for  the  predoctoral 
grant  BC050889  was  highly  successful.  The  research  conducted  yielded  exciting 
results  and  further  tested  the  first  known  reliable  model  of  ER+  breast  cancer 
metastasis  in  vivo.  This  model  and  methods  have  since  led  to  the  establishment 
of  collaborations  with  numerous  investigators  internationally.  Results  found  that 
there  are  differences  in  estrogen  responsiveness  in  breast  tumors  compared  to 
the  matched  LN  metastases,  suggesting  that  aromatase  inhibitors  and  other 
estrogen/ER  mediated  therapies  may  be  differentially  affecting  cancer  cells — 
dependent  on  their  tissue  location.  Also,  I  showed  that  laser-capture  microarray 
analysis  is  more  sensitive  and  accurate  than  whole-tissue  analysis. 
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Appendices 


1.  NIH  Biosketch. 

2.  Two  first-author  publications  from  the  past  12  months. 
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Tissue-Engineered  Bone  Is  a 
Breast  Cancer  Target 

Moreau  et  al.  _ _ _ Page  10304 

The  high  propensity  of  breast  cancer  to  metastasize  to  bone 
results  from  contributions  on  both  sides  of  the  tumor-stroma 
interface.  Identifying  the  components  in  bone  (growth  factors, 
etc.)  essential  for  metastasis  is  challenging.  To  investigate 
bone-breast  cancer  interactions,  Moreau  and  colleagues 
created  engineered  bone  by  adding  human  marrow-derived 
mesenchymal  stem  cells  to  silk  fiber  scaffolds.  Later, 
metastases  were  demonstrated  to  form  in  engineered  bone 
from  human  breast  cancer  cells  originating  at  the  orthotopic 
site  in  a  mouse  model.  Since  engineered  bone  can  be 
manipulated  in  controlled  fashion,  it  is  now  possible  to  use 
it  to  identify  bone  components  essential  for  metastasis. 

Pregnane  X  Receptor  Enhances 
Chemotherapy  in  Prostate  Cancer 

Chen  eta/.  _ Page  10361 

As  a  sensor  of  xenobiotics,  pregnane  X  receptor  (PXR) 
regulates  the  gene  expression  of  drug  metabolizing  enzymes 
and  efflux  transporters  that  can  have  profound  effects  on  the 
efficacy  of  chemotherapy.  Chen  and  colleagues  found  that 
prostate  cancer  cells  expressed  functional  PXR  which,  when 
activated,  rendered  tumor  cells  with  increased  resistance  to 
Taxol  and  vinblastine.  Ablation  of  PXR  in  prostate  cancer  cells 
enhanced  the  efficacy  of  chemotherapy.  Since  PXR  can  be 
activated  by  a  range  of  chemotherapeutics  and  herbal 
supplements,  potential  activation  of  PXR  should  be  taken  into 
consideration  when  designing  combination  chemotherapy 
for  prostate  eancer. 

Mutagenesis  Screen  Indentifies 
EGFR  Inhibitors 

Yu  et  al.  _ Page  10417 

Secondary  EGFR  mutations  commonly 
lead  to  resistance  to  the  EGFR 
inhibitors  erlotinib  and  gefitinib  in 
patients  with  non-small  cell  lung 
cancer.  Irreversible  EGFR  inhibitors, 
such  as  CL-387,785  can  overcome 
resistance  and  are  in  clinical  development.  Yu  and  colleagues 
developed  a  cell-based,  in  vitro  random  mutagenesis  screen  to 
identify  EGFR  mutations  that  confer  resistance  to  CL-387,785 


and  identified  a  number  of  novel  resistance  mutations. 

Both  an  alternative  EGFR  inhibitor  and  a  Cdk4  inhibitor 
could  overcome  such  resistance.  These  results  identify 
novel  mutations  mediating  resistance  to  irreversible  EGFR 
inhibitors  and  reveal  alternative  strategies  to  overcome  or 
prevent  resistance  in  EGFR-mutant  lung  cancers. 

Aromatase  Levels  Predict  Survival  in 
Women  with  Lung  Cancer 

Mah  etal. _ .Page  10484 

Lung  cancer  remains  the  most  lethal  type 
of  human  malignancy,  with  increased 
incidences  occurring  especially  in  women. 
Mah  and  coauthors  describe  novel  findings 
about  the  prognostic  value  of  aromatase 
expression  in  predicting  the  course  of 
non-small  eell  lung  cancer  (NSCLC).  This  group 
demonstrated  the  importance  of  estrogen  and  aromatase 
for  promoting  tumor  progression  in  a  mouse  model.  They 
further  conducted  a  population-based  study  examining 
aromatase  levels  in  human  NSCLC  using  tissue  microarray 
technology.  Lower  levels  of  aromatase  were  a  strong  predictor 
of  survival  in  women  older  than  65  years  with  earher  stage 
cancer.  The  authors  predict  that  the  use  of  aromatase 
inhibitors  may  benefit  this  patient  population. 

Lymph  Node  Microenvironment  Alters 
Sensitivity  to  Estradiol 

Harrell  et  al.  _ Page  10582 

Estrogen-responsive  breast  tumors  often  become  refractory 
to  antiestrogens  and  recur  in  organs  distant  from  the 
primary  site.  To  determine  if  the  lymph  node  (LN) 
microenvironment  alters  tumor-cell  sensitivity  to  estradiol, 
Harrell  and  colleagues  developed  a  model  to  assess  the 
hormone  responsiveness  of  xenografted  human  breast 
tumors  and  their  matched  LN  metastases.  Laser-capture 
microdissection  plus  gene  expression  profiling  confirmed 
by  immunohistochemistry,  found  that  LN  metastases  were 
estrogen  insensitive  compared  to  the  primary  tumors. 
Indeed,  some  genes  were  regulated  by  estradiol  in  the 
opposite  direction  in  tumors  and  matched  metastases. 

This  study  shows  that  the  LN  microenvironment  influences 
expression  of  therapeutic  targets  that  may  alter  cancer  cell 
survival,  recurrence,  or  response  to  treatments. 
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Abstract 

The  lymphatic  system  is  a  common  avenue  for  the  spread  of 
hreast  cancer  cells  and  dissemination  through  it  occurs  at 
least  as  frequently  as  hematogenous  metastasis.  Approximate¬ 
ly  75%  of  primary  hreast  cancers  are  estrogen  receptor  (ER) 
positive  and  the  majority  of  these  maintain  receptor  expres¬ 
sion  as  lymph  node  (LN)  metastases.  However,  it  is  unknown  if 
ER  function  is  equivalent  in  cancer  cells  growing  in  the  hreast 
and  in  the  LNs.  We  have  developed  a  model  to  assess  estrogen 
responsiveness  in  ER*  hreast  tumors  and  LN  metastases. 
Fluorescent  ER*  MCF-7  tumors  were  grown  in  ovariectomized 
nude  mice  supplemented  with  estradiol.  Once  axillary  LN 
metastasis  arose,  estradiol  was  withdrawn  (EWD),  for  1  or 
4  weeks,  or  continued,  to  assess  estradiol  responsiveness. 
On  EWD,  proliferation  rates  fell  similarly  in  tumors  and  LN 
metastases.  However,  estradiol-dependent  ER  down-regulation 
and  progesterone  receptor  induction  were  deficient  in  LN 
metastases,  indicating  that  ER-dependent  transcriptional 
function  was  altered  in  the  LN.  Cancer  cells  from  estradiol- 
treated  and  EWD  primary  tumors  and  matched  LN  metastases 
were  isolated  hy  laser  capture  microdissection.  Glohal  gene 
expression  profiling  identified  transcripts  that  were  regulated 
hy  the  tissue  microenvironment,  hy  hormones,  or  hy  hoth. 
Interestingly,  numerous  genes  that  were  estradiol  regulated 
in  tumors  lost  estradiol  sensitivity  or  were  regulated  in  the 
opposite  direction  hy  estradiol  in  LN  metastases.  We  propose 
that  the  LN  microenvironment  alters  estradiol  signaling  and 
may  contrihute  to  local  antiestrogen  resistance.  [Cancer  Res 
2007;67(21):10582-91] 

Introduction 

Interactions  between  luminal  epithelial  cells  and  the  surrounding 
microenvironment  govern  the  overall  physiology  of  the  mammary 
gland  (1),  with  development  of  epithelial  ducts  being  dependent 
on  the  activity  of  both  stromal  and  epithelial  estrogen  receptor  a 
(ERa;  ref.  2).  Integrins,  growth  factors,  and  steroid  hormone  sig¬ 
naling  pathways  all  play  an  important  part  in  maintaining  normal 
glandular  architecture  (3).  Stroma  makes  up  >80%  of  breast  volume 
and  is  composed  of  fat,  interstitial/interlobular  dense  connective 
tissue,  intralobular  loose  connective  tissue,  and  lymphatic/blood 
vessels  (4).  Luminal  epithelial  cells  associate  with  these  stromal 
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elements  as  well  as  with  basement  membrane  and  myoepithelial 
cells.  Disruption  of  this  delicate  balance  results  in  dramatic  changes 
in  both  extracellular  and  intracehular  signaling  (3). 

Seventy  percent  to  80%  of  primary  breast  tumors  express  ERa 
and  initiahy  respond  to  estradiol  withdrawal  (EWD)  or  antiestrogen 
therapies  (5).  Breast  cancer  ERa  expression  and  function  are  also 
influenced  by  the  cellular  microenvironment.  In  vitro  studies 
with  human  ER*  breast  cancer  cehs  show  that  extracellular  matrix 
proteins,  such  as  type  I  collagen  and  laminin,  can  modify  proli¬ 
ferative  responsiveness  to  estrogen  (6).  A  comparison  of  in  vitro 
and  in  vivo  models  found  that  estradiol  regulates  different  genes 
in  human  breast  tumor  xenografts  compared  with  the  identical 
cells  in  culture  (7).  These  studies  suggest  that  a  microenvironment- 
dependent  influence  on  estradiol-dependent  gene  expression  exists 
in  ER*  breast  cancers.  However,  this  issue  has  received  little  atten¬ 
tion,  and  to  the  best  of  our  knowledge,  in  vivo  experiments  that  test 
if  metastatic  microenvironments  alter  estradiol-dependent  gene 
expression  have  not  been  conducted. 

Metastasis  to  the  lymph  nodes  (LN)  is  a  key  prognostic  factor 
that  conveys  advanced  disease  status  with  the  possibility  that 
cancer  cells  have  spread  to  other  more  distant  sites.  At  diagnosis, 
30%  to  50%  of  all  breast  cancers  have  spread  to  the  draining  or 
“sentinel”  LN  (8-10),  and  if  the  primary  tumor  is  ER*,  then  ~  80% 
of  LN  metastases  retain  ER  expression  (11-13).  As  tumor  cells 
disseminate  to  LNs,  they  enter  a  distinct  microenvironment.  Here, 
they  encounter  different  supporting  elements,  such  as  reticular 
fibers,  flowing  lymph  fluid,  lymphocytes,  T  and  B  cells,  increased 
numbers  of  macrophages,  and  epithelial  reticular  cells  that  are  sur¬ 
rounded  by  adipose  tissue  (14).  It  is  unknown  if  the  differences  in 
structure,  function,  and  composition  of  the  mammary  gland  com¬ 
pared  with  the  LN  affect  ER  activities  and  estradiol-dependent  gene 
expression  in  breast  cancer  cells  that  are  located  within  these  two 
distinct  microenvironments. 

The  studies  presented  here  use  a  xenograft  model  of  fluorescent, 
ER*,  estradiol-dependent  human  breast  cancer  LN  metastasis  to 
address  this  question.  In  the  same  mouse,  LN  metastases  are  found 
to  be  transcriptionally  estradiol  insensitive  compared  with  primary 
tumors  from  which  they  originated,  despite  retention  of  ER.  We 
identify  genes  that  may  aid  in  explaining  mechanistically  why 
advanced  ER*  breast  cancers  become  nonresponsive  to  antiestro¬ 
gen  or  EWD  therapies  and  often  recur  as  metastases. 

Materials  and  Methods 

Cell  lines.  MCF-7  human  breast  cancer  cells  were  originally  purchased 
from  the  American  Type  Culture  Collection.  The  generation  of  fluorescent 
MCF-7  cells  has  been  previously  described  (15).  In  brief  ZsGreen  retroviral 
particles  (Clontech)  were  isolated  from  PT-67  packaging  cells  (Clontech), 
filtered,  and  overlaid  onto  MCF-7  cells.  Cells  were  serially  transduced  twice 
for  24  h  each  round  and  then  subjected  to  G418  selection  followed  by 
aseptic  fluorescence-associated  cell  sorting  to  isolate  a  homogeneously 
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bright  green  subpopulation.  Generation  of  PR-B-expressing  T47D  human 
breast  cancer  cells  has  been  described  previously  (16). 

Xenograft  tumor  growth  and  metastases.  All  animal  procedures  were 
done  under  a  protocol  approved  by  the  University  of  Colorado  Institutional 
Animal  Care  and  Use  Committee.  Ovariectomized  female  athymic  nu/nu 
mice  were  obtained  from  Harlan  Sprague  Dawley  at  5  to  6  weeks  of  age. 
To  establish  tumors,  animals  were  anesthetized  with  Avertin  and  injected 
into  the  opening  of  the  lactiferous  duct  of  abdominal  mammary  glands  with 
1  million  ZsGreen-expressing  MCF-7  cells  in  100  pL  of  100%  Matrigel 
(Becton  Dickinson).  Mice  were  also  implanted  with  silastic  pellets  contain¬ 
ing  171i-estradiol  (2  mg  -i-  8  mg  cellulose)  as  previously  described  (17). 
Fluorescent  whole-body  imaging  (Illumatool  9900,  Lightools  Research)  was 
done  weekly  to  determine  when  axillary  LN  metastases  arose.  For  EWD 
mice,  once  axillary  LN  metastases  were  detected,  mice  were  anesthetized 
with  Avertin  and  estradiol-releasing  pellets  were  removed  for  1  week 
some  were  removed  for  4  weeks  to  confirm  protein  expression  and  for 
proliferation  studies.  Control  mice  were  continued  on  estradiol.  Mice  were 
euthanized  by  CO2  asphyxiation  and  fluorescent  intravital  optical  imaging 
was  done  by  coupling  the  Illumatool  with  an  Olympus  SZ-61  dissecting 
microscope  and  Olympus  C-5050  digital  camera.  Fluorescent  tumors  in 
mammary  glands  and  LNs  were  removed.  Lymphatic  vessels  (LV)  draining 
fluorescent  tumor  cells  from  mammary  glands  were  also  collected. 

Serial  passaging  of  a  ZsGreen^  T47D-PR-B  LN  metastasis  was  done  by 
isolating  the  LN,  pushing  the  tissue  through  multiple  sieves,  expanding  and 
G418  selecting  human  tumor  cells  in  culture,  and  reinjecting  the  cells  into 
the  mammary  glands  of  estradiol-supplemented  recipient  mice. 

Laser  capture  microdissection  and  expression  profiling.  Mammary 
gland  tumors  (MGT)  and  LN  metastases  from  estradiol-treated  or  1-week 
EWD  mice  were  harvested,  placed  in  a  Tissue-Tek  cryomold  (EMS),  covered 
with  NegSO  frozen  section  medium  (Richard-Allan  Scientific),  and  frozen 
in  isopentane  (Sigma)  cooled  by  liquid  nitrogen.  Sections  (8  pm)  were  cut 
through  the  entire  tissue  with  a  cryotome  set  at  — 20°C.  Each  section  was 
collected  on  an  uncharged  slide  (Fisher),  and  every  10th  slide  was  stained 
for  H&E  or  processed  for  cytokeratin  18  (CK18)  immunofluorescence.  Slides 
containing  cancer  cells  were  processed  through  75%  ethanol,  water,  75% 
ethanol,  95%  ethanol,  and  100%  ethanol  for  30  s  each  and  finally  dehydrated 
with  xylene  for  1  min  before  being  immediately  subjected  to  laser  capture 
microdissection  using  an  Arcturus  Autopix.  Two  thousand  cells  for  each 
sample  were  extracted  within  30  min  of  fixation  and  frozen  at  — 70°C  until 
further  processing.  Suggested  protocols  from  Arcturus  were  followed  for 
RNA  extraction  and  two  rounds  of  amplification  using  the  PicoPure  RNA 
isolation  kit  and  RiboAmp  HS  kit  (Arcturus),  respectively.  Samples  were 
in  vitro  transcribed,  biotinylated,  hybridized  to  Affymetrix  X3P  chips  that 
interrogate  47,000  transcripts  with  61,000  probe  sets,  and  scanned  on  an 
Agilent  Bioanalyzer  2100. 

Statistical  analyses.  Raw  expression  values  for  each  tumor  and  LN  were 
downloaded  into  GeneSpring  version  7.3  (Agilent)  for  exploratory  purposes. 
Four  matched  pairs  of  tumors  and  metastatic  LNs  (2  estradiol  and  2  EWD) 
as  well  as  2  unmatched  LNs  (1  estradiol  and  1  EWD)  for  a  total  of 
10  samples  were  analyzed.  Probe  sets  that  were  absent  across  all  10  samples 
were  filtered  from  further  statistical  analysis  (18).  Remaining  probe  sets 
were  robust  multiarray  analysis  (RMA)  normalized  in  Partek  Genomics 
Suite  6.2,  and  a  three-way  mixed  model  ANOVA  analysis  was  done  to  iden¬ 
tify  differentially  expressed  genes.  Four  groups  were  compared:  estradiol- 
treated  MGT,  estradiol-treated  LN,  EWD  MGT,  and  EWD  LN.  These  four 
groups  were  divided  based  on  the  two  levels  of  hormone  treatment 
(estradiol  or  EWD)  and  the  two  levels  of  tissue  type  (MGT  or  LN).  Hormone 
and  tissue  were  therefore  defined  as  main  factors  in  the  model.  Interac¬ 
tion  effects  among  hormone  and  tissue  were  also  included.  The  random 
factor  animal  (which  is  nested  in  hormone)  was  added  to  account  for  the 
fact  that  multiple  tissues  (MGT  and  LN)  were  taken  from  the  same  animal. 
In  total,  1,570  transcripts  were  identified  as  having  significant  differences 
among  the  four  group  means  (.Pmodei  <  0.05)  and  P  <  0.05  for  main  or 
interaction  effects.  To  further  identify  differences  in  hormonal  regulation 
among  the  1,570  differentially  expressed  genes,  they  were  subjected  to  linear 
contrasts  to  make  two  specific  comparisons:  EWD  versus  estradiol-treated 
MGT  and  EWD  versus  estradiol-treated  LN.  Bonferroni’s  correction  for 


multiple  tests  was  then  applied  to  these  two  comparisons  and  the  level  of 
significance  was  changed  to  0.025. 

Gene  networking.  To  define  relationships  among  genes  differentially 
regulated  by  estradiol  in  MGTs  and  LN  metastases,  pathway  analyses  were 
done  using  Ingenuity  software.  Probe  set  identification  numbers  were  used 
for  gene  identification. 

Immunohistochemistry.  Organs  were  fixed  in  4%  paraformaldehyde 
overnight,  paraffin  embedded,  and  cut  into  4-pm  sections.  After  high- 
temperature  antigen  retrieval  in  citrate  buffer,  immunohistochemistries 
were  done  on  estradiol-treated  tissues  and  1-week  or  4-week  EWD  mice 
using  primary  antibodies  directed  against  ER  (1:100,  SPl;  Lab  Vision/ 
NeoMarkers),  progesterone  receptor  (PR;  1:500,  1294;  Dako),  CK18  (1:200; 
Calbiochem),  CD44  (1:200;  Ab4;  Lab  Vision/NeoMarkers),  caveolin  1  (CAVl; 
1:125;  Epitomics),  and  cathepsin  D  (CTSD;  1:125;  Epitomics)  for  1  h  at 
room  temperature.  Bound  primary  antibodies  were  detected  with 
horseradish  peroxidase  (HRP)-conjugated  goat  anti-mouse  and/or  goat 
anti-rabbit  secondary  antibodies  (Envision  HRP;  Dako)  reacted  with  3,3'- 
diaminobenzidine  (Dako).  Sections  were  counterstained  with  hematoxylin 
and  mounted  with  Permount  (Fisher).  For  immunofluorescence  microscopy, 
goat  anti-rabbit  Alexa  Fluor  555  (red)  secondary  antibodies  were  used 
(1:200,  Alexa  Fluors;  Invitrogen)  and  sections  were  counterstained  with 
4',6-diamidino-2-phenylindole  (DAPI)  in  methanol  and  mounted  with  Gel/ 
Mount  (Biomeda). 

Results 

EWD  model.  MCF-7+ZsGreen  tumors  require  estradiol  for 
growth  and  LN  metastasis  (15).  To  generate  LN  metastases, 
bilateral  tumors  were  grown  in  the  abdominal  mammary  glands 
of  ovariectomized  immunocompromised  mice  that  were  supple¬ 
mented  with  an  estradiol-releasing  pellet  (Fig.  1.4).  Weekly  tumor 
measurements  were  recorded  with  a  digital  caliper.  Spread  of 
cancer  cells  to  axillary  LNs  was  determined  by  weekly  fluorescent 
whole-body  imaging  (Fig.  IB).  To  define  estradiol-regulated  genes 
(7),  once  LN  metastases  were  detected,  the  estradiol-releasing 
pellet  was  withdrawn  for  1  week  in  half  of  the  mice  (EWD)  or  left  in 
place  in  the  remaining  mice.  To  assess  proliferation  rate,  EWD  was 
extended  for  4  weeks  in  tumors  and  LN  metastases.  Bromodeox- 
yuridine  labeling  and  mitotic  indices  indicated  a  similar  decrease 
of  ~  50%  at  both  sites.  At  necropsy,  assessment  of  the  uterine 
mass  confirmed  the  presence  or  absence  of  systemic  estradiol,  with 
1-  and  4- week  EWD  uteri  exhibiting  significant  decreases  in  mass 
compared  with  estradiol-treated  uteri.  To  locate  cancer  cells  within 
LNs,  the  entire  LN  was  sectioned  and  every  10th  slide  was  pro¬ 
cessed  for  CK18  immunohistochemistry  (Fig.  1C).  The  CK18  map 
and  the  unique  morphology  of  MCF-7  cells  compared  with 
lymphocytes  (Fig.  ID)  allowed  for  identification  of  the  cancer  cells 
in  serial  sections.  We  then  marked  and  laser  captured  pure  MCF-7 
cell  populations  from  matched  tumors  and  LN  metastases. 

Estrogen  insensitivity  in  LN  metastases:  ER  and  PR 
expression.  Estradiol  treatment  leads  to  down-regulation  of  ER 
protein  levels  in  tumors  (19,  20).  To  characterize  these  effects  in 
estradiol-free,  estradiol-treated,  and  EWD  conditions,  the  MCF-7 
tumors,  LV  emboli,  and  LN  metastases  were  probed  for  ER  and 
PR  by  immunohistochemistry.  In  Fig.  2,  the  tissues  from  each 
treatment  group  were  from  the  same  mouse.  Highest  ER  levels 
were  present  in  tumors  that  did  not  receive  supplemental  estradiol 
(— E).  Seventy-three  percent  of  MCF-7  cells  expressed  high  levels 
of  ER  and  no  cells  contained  PR  (Fig.  2,  top).  In  the  presence 
of  estradiol  (+E),  ERs  were  down-regulated  by  >40%  in  the  tumor 
(Fig.  2,  middle)  but  decreased  by  <20%  in  LV  tumor  emboli  and  LN 
metastasis  of  the  same  mouse.  Ligand-dependent  down-regulation 
of  ER  is  required  for  its  transcriptional  activity  (20),  including  PR 
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Figure  1.  MCF-7+ZsGreen  xenograft  tumors  and  LN  metastases. 
MCF-7+ZsGreen  tumors  were  established  in  mammary  glands  of  ovariectomized 
nu/nu  mice  supplemented  with  an  estradiol  pellet,  and  LN  metastases 
were  allowed  to  develop.  A,  fluorescent  whole-body  image  of  bilateral 
MCF-7-rZsGreen  tumors.  B,  fluorescent  whole-body  image  of  MCF-7-rZsGreen 
metastasis  to  the  left  axillary  LN.  C,  fluorescence  immunohistochemistry  for 
CK18  (red)  and  DAPI  (blue)  of  an  axillary  LN  metastasis.  Magnification,  x4. 

D,  white  light  laser  capture  microscope  image  of  an  axillary  LN  metastasis 
showing  MCF-7  cells  or  lymphocytes  (LN).  Magnification,  x40. 


induction.  Indeed,  PRs  were  expressed  in  -i-E  tumors  but  were  only 
rarely  found  in  -i-E  LV  emboli  (Fig.  2,  arrows)  or  in  LN  metastases. 
On  EWD,  ERs  were  restored  to  control  (— E)  levels  in  aU  tissues,  and 
PRs  were  absent  (Eig.  2,  bottom). 


Microenvironment  versns  clonal  selection.  To  determine 
whether  the  increased  ER  expression  observed  in  LN  metastases 
was  due  to  clonal  selection  of  ER*  cells  from  the  primary  tumor,  or 
to  an  effect  on  ER  at  the  metastatic  microenvironment,  tumor  cells 
were  isolated  from  the  axillary  LN  metastasis  of  a  donor  mouse, 
expanded  in  culture,  and  reimplanted  in  the  mammary  gland  of  an 
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Figure  2.  Estradiol  regulation  of  ER  and  PR  levels.  ER  and  PR  expression 
In  MCF-7  cells  from  matched  tumors,  LVs,  and  LNs.  MCF-7-rZsGreen  tumors 
were  established  in  ovariectomized  mice  supplemented  with  a  placebo  pellet 
(-E),  an  estradiol-releasing  pellet  (+E),  or  an  estradiol-releasing  pellet  until  LN 
metastases  arose,  after  which  the  estradiol  pellet  was  withdrawn  for  1  wk 
(EWD).  The  primary  tumors,  efferent  LVs,  and  draining  LNs  from  the  same 
mouse  were  dissected  and  prepared  for  histology.  Paraffin-embedded  sections 
were  probed  with  rabbit  anti-ER  or  mouse  anti-PR  primary  antibodies  coupled 
with  goat  anti-rabbit  (ER)  or  goat  anti-mouse  (PR)  HRP-conjugated  secondary 
antibodies.  Background  staining  of  erythrocytes  (e).  Arrows  point  to  the  LV. 
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Figure  3.  Aberrant  ER  regulation  in  tumor 
cells  located  in  the  LN  microenvironment 
and  global  gene  expression  profiling 
of  tumors  and  matched  LN  metastases. 

A,  ER  immunohistochemistry  of  estradiol- 
freated  T47D  tumors  and  their  matched  LN 
metastases.  An  estradiol-treated  MGT  in 
mouse  1  developed  LN  metastases. 

Cells  from  the  LN  metastasis  of  mouse  1 
were  expanded  in  culfure  and  reinjected 
into  the  mammary  gland  of  a  second 
estradiol-treafed  recipienf  mouse  2, 
which  also  yielded  LN  mefasfases.  MGTs 
and  LN  metastases  from  both  mice 
were  paraffin  embedded,  secfioned,  and 
stained  for  ER  by  immunohistochemistry. 

B,  unsupervised  hierarchical  clustering 
from  confinuously  estradiol-freafed  (+E) 
or  1-wk  EWD  tumors  and  matched 

LN  metastases.  Red,  up-regulated 
transcripts;  blue,  down-regulated 
transcripts;  gray,  relatively  unchanged 
transcripts.  C,  ANOVA  analysis  identified 
1 ,570  genes  with  tissue  effects,  hormone 
effects,  and/or  interaction  effects. 
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estradiol-treated  recipient  mouse.  This  tumor  also  yielded  a  LN 
metastasis.  Figure  3A  (top)  shows  ER  immunohistochemistry  of 
the  estradiol-treated  primary  MGT  and  its  matched  LN  metastasis 
from  the  donor  mouse  (mouse  1).  ERs  in  the  tumor  grown  in  the 
mammary  gland  and  LN  of  the  recipient  mouse  (mouse  2)  are 
shown  in  Fig.  3A  (bottom).  There  was  a  2-fold  increase  in  ER 
content  in  the  cancer  cells  populating  the  LN  of  mouse  1  and  a 
similar  2-fold  increase  in  the  LN  ER  of  mouse  2.  These  data  suggest 
that  the  LN  microenvironment  rather  than  clonal  selection 
modifies  ER  expression  levels  in  response  to  estradiol. 

Global  gene  expression  profiling.  To  determine  if  the  apparent 
estradiol  insensitivity  of  LN  metastases  as  assessed  by  ER  down- 
regulation  and  failure  of  PR  induction  extends  to  other  estradiol- 
dependent  genes,  global  gene  expression  profiling  was  done.  To 
ensure  that  only  human  tumor  cells  were  evaluated,  pure  cell  popu¬ 
lations  were  isolated  by  laser  capture  microdissection.  Estradiol- 
treated  and  1-week  EWD  tumors  and  their  matched  LN  metastases 
were  analyzed  to  define  the  subset  of  estradiol-regulated  genes  (7). 
Unsupervised  hierarchical  clustering  showed  that  EWD  had  a  lesser 
effect  on  genes  in  LN  metastases  than  on  genes  in  MGTs,  with  the 
estradiol-treated  and  EWD  LN  metastases  clustering  closest 


together  (Fig.  3B).  Interestingly,  the  EWD  tumors  clustered  closer 
to  the  LNs  than  to  the  estradiol-treated  tumors. 

To  identify  genes  with  differential  expression  patterns,  a  three- 
way  mixed  model  ANOVA  was  done  for  genes  that  were  present 
in  at  least  one  condition.  Genes  with  the  greatest  differences  in 
expression  due  to  the  hormone  treatment,  tissue  type,  interaction 
between  hormone  and  tissue,  or  within  and  between  animals 
(^modei  <  0-05)  and  with  significant  (P  <  0.05)  main  (tissue  and/ 
or  hormone  effects)  or  interaction  effects  were  sorted  to  yield  a 
total  of  1,570  differentially  expressed  transcripts.  Significant  main 
effects  indicate  simple  differences  between  the  two  hormone  and/ 
or  two  tissue  groups,  an  example  being  that  all  EWD  tissues  have 
higher  gene  expression  levels  than  all  estradiol-treated  tissues 
(regardless  of  tissue  type)  or  that  all  LN  tissues  have  higher  gene 
expression  levels  than  all  MGT  (regardless  of  hormone  treatment). 
Significant  interaction  effects,  on  the  other  hand,  indicate  com¬ 
plex  differences  in  which  the  results  of  one  factor  differ  depending 
on  the  levels  of  the  other  factor.  For  example,  a  gene  that  was 
expressed  at  higher  levels  in  estradiol-treated  LNs  than  MGTs  but 
in  EWD  tissues  had  the  opposite  pattern  with  higher  levels  in 
MGTs  than  LNs. 
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Figure  4.  Genes  that  are  tissue  and 
hormone  regulated  with  or  without 
interaction  effects.  RMA  normalized  mean 
RNA  expression  values  with  SE.  Solid 
lines,  data  from  estradiol-treated  tissues; 
dashed  lines,  data  from  EWD  tissues. 

A,  tissue-regulated  genes.  B,  hormone- 
regulated  genes.  MT1F,  metailothionein 
1 F.  C,  genes  that  are  tissue  and  hormone 
regulated.  D,  genes  with  interaction  effects 
that  are  differentially  regulated  by  estradiol 
in  the  two  tissues  include  CAV1,  MUC5AC, 
TP53,  BCL2,  CDKN1A,  and  CD44. 


Of  the  1,570  differentially  expressed  transcripts,  more  had 
increased  average  expression  in  LN  metastases  (70%)  than  in 
primary  tumors  (30%).  In  addition,  on  average,  genes  were  more 
highly  expressed  in  EWD  tissues  (64%)  than  in  estradiol-treated 
ones  (36%),  suggesting  that  estradiol  affects  approximately  two 
thirds  of  genes  through  suppression  of  RNA  levels.  Among  the 
1,570  transcripts,  386  had  only  significant  [P  <  0.05;  Fig.  3C)  tissue 
effects,  44  had  only  significant  hormone  effects,  and  91  were 
affected  by  both  the  tissue  microenvironment  and  the  presence  of 
hormones  (Fig.  3C).  However,  the  majority  of  genes,  1,049  of  1,570, 
exhibited  a  significant  interaction  effect.  The  expression  values 
of  these  genes  could  not  be  simply  described  as  due  to  either  of 
the  main  effects  but  rather  due  to  an  interrelationship  between 
them.  Most  of  these  genes  may  also  have  significant  main  effects  as 
shown  by  intersection  in  the  Venn  diagram,  but  interaction  effects 
take  precedence  over  these  main  effects,  and  the  main  effects 
cannot  be  correctly  interpreted  without  acknowledging  the  inter¬ 
action  effects  (21).  For  these  genes,  expression  levels  in  the  two 
tissues  (MGT  and  LN)  were  influenced  differently  by  the  type  of 


hormone  treatment  (estradiol  and  EWD).  A  subset  (145  of  1,049)  of 
these  genes  exhibited  only  interaction  effects. 

Main  effects  and  interaction  effects  can  be  visualized  by 
graphical  representation  with  levels  of  one  factor  on  the  X  axis 
and  separate  lines  representing  levels  of  the  other  factor,  as  in 
Fig.  4.  Lines  that  are  parallel  to  one  another  with  steep  slopes  or 
separation  indicate  main  effects,  whereas  nonparaUel  lines  indicate 
interaction  effects.  Figure  4A  shows  two  examples  of  genes  with 
significant  tissue  effects.  Cortactin  transcripts  had  significantly 
higher  expression  levels  in  LN  metastases,  regardless  of  hormone, 
as  did  the  majority  (290  of  386)  of  the  tissue-regulated  genes  (other 
examples  include  MAPK  and  cyclin  F).  CXCL14  is  one  of  the  96  of 
386  tissue-regulated  genes  whose  levels  were  higher  in  the  primary 
tumors  compared  with  LN  metastases. 

Figure  45  shows  two  examples  of  genes  with  a  significant  hor¬ 
mone  effect.  The  expression  of  cathepsin-D  (CTSD)  transcripts  was 
dependent  on  estradiol  regardless  of  the  microenvironment.  Other 
genes  with  similar  expression  patterns  were  insulin-like  growth 
factor-I  receptor  {IGF-IR),  cadherin-15,  and  defensin  p-1.  Transcript 
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levels  of  metallothionein  IF  were  suppressed  by  estradiol  regardless 
of  the  microenvironment.  Other  examples  of  estradiol-suppressed 
genes  are  cathepsin  F  and  BLNK. 

A  subset  of  differentially  expressed  genes  (91  of  1,570)  had  RNA 
expression  levels  that  were  both  hormone  and  tissue  regulated  as 
shown  in  Fig.  4C.  ERBB3  was  significantly  suppressed  by  estradiol 
regardless  of  the  microenvironment,  but  its  expression  was  also 
increased  in  LN  metastases  compared  with  MGTs  of  both  hormone 
treatments.  SERPl  was  also  significantly  increased  in  LN 
metastases  compared  with  MGTs  regardless  of  the  presence  of 
estradiol  and  was  also  significantly  up-regulated  by  estradiol  in 
both  tissues. 

The  majority  of  differentially  expressed  genes  (1,049  of  1,570)  had 
a  significant  interaction  effect,  exhibiting  an  interrelationship 
between  tissue  and  hormonal  factors.  Figure  4D  shows  examples  of 
a  few  such  genes.  Transcripts  for  BCL2,  CAVl,  CD44,  CDKNIA, 
MUC5AC,  and  tumor  promoter  53  (TP53)  all  exhibit  interaction 
effects  and  are  affected  differently  by  estradiol  in  tumors  and  LN 
metastases.  MUC5AC  and  TP53  are  examples  of  genes  exhibiting 
main  effects  (hormone  and  tissue)  but  with  the  interaction  term 
predominating.  Main  effects  of  CDKNIA  {p21)  cancel  each  other 
out,  leaving  the  dominant  role  of  interaction  effects  evident. 

Immunohistochemistry.  Several  of  the  genes  that  were 
identified  by  the  expression  profiling  analysis  were  investigated 
further  at  the  protein  level  by  immunohistochemistry  (Fig.  5).  The 
estradiol-regulated  gene  CTSD  was  expressed  in  80%  to  90%  of 
estradiol-treated  MCF-7  primary  tumor  cells,  and  EWD  reduced 
this  number  to  50%  to  60%.  In  the  LN  metastases  in  the  presence  of 
estradiol,  50%  to  60%  of  cancer  cells  expressed  CTSD;  this  fell  to 
10%  to  15%  on  EWD  (Fig.  5,  top).  CAVl  protein  was  unaffected  by 
EWD  but  was  expressed  only  in  MGTs  and  not  in  LN  metastases 
(Fig.  5,  middle).  In  estradiol- treated  tumors,  CD44  was  heteroge¬ 
neously  expressed;  on  EWD,  nonsignificant  increases  in  the 
percentage  of  CD44*  cells  were  recorded.  Interestingly,  the  opposite 
pattern  was  observed  in  LN  metastases.  The  extensive  over¬ 
expression  of  CD44  in  estradiol-treated  cancer  cells  metastatic  to 
LNs  (15)  is  dramatically  reduced  by  EWD  (Fig.  5,  bottom). 

Estrogen-regulated  genes.  Two  linear  contrasts  were  done 
among  the  1,570  differentially  expressed  transcripts  comparing 
estradiol-treated  with  EWD  tumors  and  estradiol-treated  with 
EWD  LN  metastases  (Fig.  6A  and  B).  This  identified  transcripts  that 
were  estradiol  regulated  in  MGTs  compared  with  LN  metastases. 
There  were  273  significantly  hormonally  up-regulated  or  down- 
regulated  genes  (fold  change  >1.3;  P  <  0.025;  Fig.  6A).  Of  these,  256 
were  estradiol  dependent  in  tumors,  41  were  estradiol  dependent  in 
LN  metastases,  and  24  maintained  estradiol  sensitivity  in  both 
MGTs  and  LN  metastases  (Fig.  6A).  These  data  confirm  the  con¬ 
clusions  drawn  from  Figs.  2  and  SA  that  tumor  cells  in  the  LN 
microenvironment  are  relatively  estradiol  resistant. 

To  identify  genes  with  any  absolute  difference  in  expression  due 
to  estradiol  in  MGTs  and  LN  metastases,  the  fold  change  criterion 
was  removed  from  the  two  linear  contrasts.  This  increased  the 
total  list  of  signihcantly  {P  <  0.025)  hormone-regulated  genes  horn 
273  to  1,176  (Fig.  6B).  These  additional  903  genes  were  up-regulated 
or  down-regulated  with  fold  changes  >1.0  and  <1.3,  and  in  total, 
there  were  707  estradiol-regulated  genes  in  MGTs,  177  estradiol- 
regulated  in  LN  metastases,  and  292  estradiol-regulated  in  both 
sites.  Seventy-five  of  these  transcripts  were  regulated  by  estradiol 
in  opposite  directions  in  MGTs  and  LN  metastases  (Fig.  6C). 
BCL2,  CD44,  and  CDKNIA  are  examples  of  such  genes  (Fig.  AD). 
CDKNIA,  also  known  as  p21/WAF,  encodes  a  potent  cyclin- 
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Figure  5.  Immunohistochemical  analyses  of  proteins  that  are  hormone  and 
tissue  regulated.  Estradiol-treated  and  EWD  tumors  and  matched  LN 
metastases  from  the  same  mouse  were  paraffin  embedded,  sectioned,  and 
probed  with  CTSD,  CAVl,  or  CD44  primary  antibodies  coupled  with  goat 
anti-rabbit  or  goat  anti-mouse  HRP-conjugated  secondary  antibodies.  Top, 
estradiol-regulated  CTSD;  middle,  CAVl  protein  regulation  by  the  tissue 
microenvironment;  bottom,  differential  estradiol  regulation  of  CD44  expression 
in  tumors  and  LN  metastases. 


dependent  kinase  inhibitor  and  regulator  of  cell  cycle  progression 
at  Gi.  Its  expression  is  tightly  controlled  by  TP53,  one  of  our  set 
of  microenvironment-regulated  genes  (Fig.  AD).  CD44  (Fig.  5)  is  a 
putative  breast  cancer  stem  cell  marker  (22)  and  may  target  tumor 
cells  to  LNs  (15). 

Network  analysis.  Ingenuity  data  analysis  was  conducted  to 
determine  if  the  genes  that  had  interaction  effects  and  were 
differentially  regulated  by  estradiol  in  MGTs  and  LN  metastases  are 
also  involved  in  important  functional  networks.  Twenty-seven 
unique  gene  networks  that  contained  at  least  one  differentially 
regulated  gene  were  identified.  Supplementary  Fig.  SI  shows  a 
network  of  33  genes,  100%  of  which  are  differentially  estradiol 
regulated  in  MGTs  and  LN  metastases.  This  network  has  TP53  at 
its  center  and  functions  in  cellular  responses  to  “therapeutics, 
cell  death,  and  connective  tissue  disorders.”  Most  of  these  genes 
(25  of  33)  were  estradiol  regulated  in  MGTs  but  lost  estradiol  sen¬ 
sitivity  in  LN  metastases  (orange  outline).  Two  genes  were  estradiol 
regulated  in  LN  metastases  but  not  in  tumors  (red  outline),  and  the 
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remaining  seven  genes  were  regulated  in  the  opposite  direction  by 
estradiol  in  MGTs  and  LN  metastases  (blue  outline). 

Discussion 

Gene  expression  changes  in  MGTs  and  LN  metastases.  Two 

major  theories  have  been  advanced  to  explain  successful  cancer 
metastasis,  (a)  The  “rare  tumor  cell  hypothesis”  states  that  only 
certain  cells  within  tumors  are  capable  of  migrating  and/or 
growing  at  distant  sites.  These  rare  cells  may  be  stem  cells,  such  as 
those  that  express  CD4f  (22).  On  the  other  hand,  they  may 
represent  rare  cells  within  a  population  of  genetically  diverse  and 
unstable  cells  that  originated  with  or  acquired  the  appropriate 
genetic  information  needed  to  seed  at  specific  metastatic  sites  (23). 
[b)  An  alternative  hypothesis  posits  that  metastatic  capacity  is  an 
inherent  rather  than  an  acquired  feature  of  breast  cancers  (24,  25) 
and  is  based  on  data  suggesting  that  metastatic  capacity  is 
encoded  early  in  the  majority  of  cells  within  a  tumor.  In  support 
of  this  idea  are  studies  showing  that  tumors  and  their  LN  meta¬ 
stases  retain  similar  genetic  signatures,  including  aggressiveness 
and  proliferation  markers  (24,  26,  27).  One  study,  for  example, 
found  only  27  changed  genes  between  lung  tumors  and  their  LN 
metastases  (28).  Weigelt  et  al.  (27)  identified  a  70-gene  cluster  in 
primary  breast  tumors  that  predicts  likelihood  of  subsequent 
metastases  because  the  unique  expression  profiles  of  these  genes 
were  maintained  in  LNs  and  more  distant  metastases.  Another 
study  found  no  “shared  genes”  between  primary  tumors  that  were 
predictive  of  LN  metastases  but,  anywhere  between  3  and  149 
genes,  were  “antiexpressed”  or  expressed  in  opposite  directions  in 
tumors  versus  LN  metastases  (29). 

Our  data  provide  evidence  in  support  of  both  theories.  In  these 
studies,  we  analyze  the  changes  in  gene  expression  induced  by 
estradiol  as  revealed  by  estradiol  removal.  Fewer  than  5%  of  all 
genes  analyzed  were  differentially  expressed  in  cancer  cells  within 
the  mammary  gland  compared  with  the  LN.  Although  this  might  be 
considered  a  subtle  difference,  it  nevertheless  represents  differen¬ 
tial  regulation  of  >1,000  genes  in  these  two  microenvironments. 
These  genetic  differences  could  be  sufficient  to  define  a  “metastatic 
stem  cell”  (30),  expressing  the  adaptive  proteins  required  for 
successful  colonization  and  growth  of  a  subpopulation  in  a  meta¬ 
static  microenvironment,  as  envisioned  by  the  rare  tumor  cell 
hypothesis.  Elevated  expression  levels  of  the  hyaluronan  receptor 
CD44,  a  putative  cancer  stem  cell  marker,  in  estradiol-treated  LN 
metastases  suggest  such  selectivity  in  our  models.  On  the  other 
hand,  our  observations  that  cancer  cell  emboli  move  in  bulk 
through  LVs  to  LNs  (15),  coupled  with  our  findings  here  and  those 
of  others  (24,  27,  29)  of  modest  genetic  changes  between  primary 
tumors  and  LN  metastases,  could  be  interpreted  as  supporting 
the  view  that  the  majority  of  primary  tumor  cells  arise  either  with 
or  without  the  predisposition  for  LN  spread.  If  true,  the  gene 
expression  changes  we  observe  could  be  attributed  to  the  plasticity 
that  all  tumor  cells  exhibit  when  growing  in  two  distinct  micro¬ 
environments.  Evidence  supporting  these  microenvironmental 
effects  on  ER  expression  and  function  was  found  when  a  down- 
regulation  “resistant”  ER“  LN  metastasis  regained  ER  down- 
regulation  capacity  when  grown  in  the  mammary  gland  of  a 
recipient  mouse  (Fig.  3A). 

ER  function  in  metastases.  Beatson  (31)  reported  >100  years 
ago  that  an  essential  feature  of  breast  cancer  is  continuous  and 
excessive  growth  of  the  epithelium,  which  invades  the  surround¬ 
ing  tissues,  spreads  along  the  LVs,  passes  from  one  set  of  LNs  to 
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another,  and  eventually  forms  deposits  in  distant  organs,  steps  that 
occur  more  rapidly  and  prove  to  be  more  quickly  fatal  in  younger 
patients.  Knowing  that  the  ovaries  produced  factors  required  for 
lactation,  and  believing  “that  all  pathological  changes  are  merely 
modified  physiological  ones,”  Beatson  was  the  first  to  find  that 
removal  of  the  ovaries  and  the  estradiol  they  elaborate  reduced 
breast  tumor  size  and  extended  the  survival  of  young  women  with 
advanced  breast  cancer. 

Seventy  percent  to  80%  of  primary  breast  cancers  are  ER*,  and 
~80%  of  LN  metastases  retain  the  receptors  (11-13,  32).  If  such 
tumors  are  treated  with  antiestrogens  or  EWD  therapies,  20%  to 
50%  acquire  resistance  (33),  often  recurring  as  LN  metastases. 
Premenopausal  women  with  LN  metastases  are  less  responsive  to 
tamoxifen  treatment  than  LN-negative  patients  (34),  and  nodal 
involvement  is  a  significant  predictor  of  local  relapse  (35)  with  a 
5.6-fold  higher  relative  risk  of  local  recurrence  than  women  with 
clear  LNs.  Preoperative  chemotherapy  causes  significantly  less 
damage  to  LN  metastases  in  comparison  with  primary  tumors  (36). 
Based  on  these  data,  we  hypothesized  that  ER  effects  on  gene 
expression  would  be  altered  by  the  LN  microenvironment.  To  the 
best  of  our  knowledge,  no  previous  studies  have  defined  estradiol- 
regulated  genes  in  tumors  compared  with  their  matched  LN  meta¬ 
stases.  Identifying  these  genes  is  an  important  step  toward 
elucidating  mechanisms  that  contribute  to  breast  cancer  progres¬ 
sion  and  estradiol  resistance  in  advanced  disease. 

In  the  face  of  continuous  estradiol  treatment  in  two  separate  ER* 
breast  cancer  cell  lines,  higher  levels  of  ER  were  observed  in  LN 
metastases  than  in  matched  MGTs  (Eigs.  2  and  3A).  The  partial 
failure  of  ligand-dependent  ER  down-regulation  in  LN  metastases 
suggests  that  the  receptors  are  relatively  estradiol  insensitive  in  the 
LN  microenvironment.  Reduced  PR  expression  in  estradiol-treated 
LN  metastases  supports  this  conclusion,  as  does  global  gene 
expression  profiling,  which  showed  that  estradiol  had  less  influence 
on  gene  expression  in  LN  metastases  than  in  matched  primary 
MGTs  (Pigs.  3  and  6).  We  have  now  identified  numerous  genes  that 
are  estradiol  regulated  in  primary  tumors  but  not  in  LN  metastases. 
Examples  include  MUC5A  and  TP53,  which  were  estradiol  sensitive 
in  MGTs  but  lost  this  regulation  in  the  LNs.  Additionally,  75  genes 
were  regulated  in  the  opposite  direction  by  estradiol  in  MGTs 
compared  with  their  LN  metastases  (Pig.  6C),  21  of  which, 
including  BCL2  and  CDKNIA  (Fig.  4),  were  previously  reported  to 
be  estradiol  regulated  (7,  37,  38).  CD44  (Figs.  4  and  5)  was  also 
regulated  in  opposite  directions  at  the  two  sites.  It  is  one  of  two  cell 
surface  markers  that  differentiate  tumorigenic  from  nontumori- 
genic  breast  cancer  cells  (22).  This  putative  stem  cell  and  aggres¬ 
siveness  marker  (39,  40)  is  an  example  of  how  EWD  differentially 
alters  the  expression  of  a  potential  therapeutic  target  in  primary 
tumors  compared  with  LN  metastases. 

Mechanisms  of  estradiol  resistance  in  LNs.  Possible  explan¬ 
ations  for  the  differential  estradiol  sensitivity  of  MGTs  versus  LN 
metastases  include  increased  availability  of  growth  factors  in  LNs 
that  directly  interfere  with  ER  transcription  (41),  diminished 
ubiquitination  and  loss  of  proteasome-mediated  protein  down- 
regulation  in  LNs  (42,  43),  loss  of  CAVl  and  subsequent  membrane 
signaling  (44),  and  reduced  estradiol  availability  to  LNs.  With 
regard  to  growth  factors,  IGFs  and  epidermal  growth  factors  (EGF) 
activate  the  phosphatidylinositol  3-kinase-AKT-mammalian  target 
of  rapamycin  pathways,  which  alter  Jun  and  Eos  activities  at  the 
PR  promoter  and  decrease  its  ER-mediated  transcription  (45). 
LN  stromal  cells  produce  both  IGF  and  EGF  (41),  which  may 
explain  the  reduced  PR  induction  as  well  as  the  reduced  estradiol 
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Figure  6.  Pairwise  comparisons  to  identify  estradiol-reguiated  transcripts  in  tumors  (7),  LN  metastases,  or  both.  A,  transcripts  reguiated  in  a  significant  manner 
(P  <  0.025)  by  estradiol  >1 .3-foid.  B,  transcripts  reguiated  in  a  significant  manner  (P  <  0.025)  by  estradioi  regardiess  of  foid  ieveis.  C,  iist  of  75  transcripts  that  are 
reguiated  in  a  significant  manner  in  opposite  directions  by  estradiol  in  MGTs  compared  with  LN  metastases.  /,  increases;  D,  decreases  with  estrogen. 


regulation  of  other  target  genes  at  this  site.  Increased  cortactin 
levels  (Fig.  4)  inhibit  ubiquitin-mediated  degradation  of  EGF 
receptors  following  their  internalization,  resulting  in  sustained 
ligand-induced  EGF  signaling  (42,  43).  It  is  possible  that  similar 


ubiquitin-dependent  mechanisms  explain  the  ER  down-regulation 
defects  we  observe  in  LNs  (Fig.  2;  ref.  15).  Supporting  this  hypo¬ 
thesis,  nine  transcripts  known  to  be  involved  in  the  ubiquitination 
process  were  differentially  expressed  in  MGTs  and  LN  metastases. 
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This  includes  two  known  enzymes  that  function  in  the  ubiquitina- 
tion  process,  UBE2D3  and  USP7,  which  were  increased  in  LN 
metastases  and  also  lost  estradiol  sensitivity  in  LN  metastases 
(Supplementary  Fig.  SI).  CAVl,  a  scaffold  protein  that  transduces 
ER  signaling  at  the  plasma  membrane  (44),  was  decreased  in  LNs, 
possibly  explaining  loss  of  ER  responsiveness.  Finally,  although  it  is 
possible  that  the  diminished  estradiol  sensitivity  of  LN  meta¬ 
stases  is  due  to  an  overall  decrease  of  estradiol  levels  within  LNs, 
this  is  highly  unlikely  given  the  number  of  estradiol-regulated  genes 
whose  levels  are  maintained  or  increased  in  LN  metastases  (Fig.  4). 

Networks.  We  speculated  that  genes  differentially  regulated  in 
MGTs  versus  LN  metastases  may  contribute  to  tumor  aggressive¬ 
ness,  recurrence,  and/or  growth  in  LNs.  We  therefore  analyzed  the 
functional  networks  characterizing  these  genes.  Of  interest  was  a 
33-gene  network  composed  entirely  of  transcripts  regulated 
differently  by  estradiol  in  MGTs  and  LNs  (Supplementary  Fig.  SI), 
which  included  CDNKIA,  TP53,  and  BCL2-linked  signaling  path¬ 
ways.  In  this  network,  the  majority  of  the  genes  were  significantly 
regulated  by  estradiol  in  tumors  while  losing  estradiol  sensitivity  in 
LN  metastases  or  suppressed  by  estradiol  in  tumors  but  increased 
by  estradiol  in  LN  metastases.  These  findings  are  interesting  given 
the  fact  that  70%  of  the  1,570  global  genes  differentially  expressed 
in  the  two  sites  were  increased  in  LN  metastases  compared  with 
MGTs,  whereas  64%  of  the  differentially  estradiol-regulated  genes 
were  suppressed  by  estradiol.  This  corroborates  previous  data  with 
MCF-7  cells  showing  that  70%  of  genes  are  down-regulated  by 
estradiol  treatment  (38)  and  is  consistent  with  an  overall  lifting  of 
the  suppressive  effects  of  estradiol  in  MCF-7  tumor  cells  growing  in 
LNs.  It  is  possible  that  these  results  do  not  completely  mimic 
effects  in  a  patient  with  an  intact  immune  system.  The  LNs  of  nude 
mice  lack  or  have  reduced  numbers  of  T  cells  (46).  However,  they 
have  B  cells,  natural  killer  cells,  reticular  cells,  and  macrophages  as 
well  as  a  nodal  stromal  network  and  architecture.  Because  different 


substrata  are  known  to  alter  growth  and  function  of  ER*  human 
breast  cancer  cells  in  vitro  (6),  it  is  likely  that  the  LN  micro¬ 
environment,  including  its  unique  architecture,  cytokines,  and 
growth  factors,  modifies  ER  function.  Given  the  experimental 
nature  of  our  models,  these  findings  stress  the  need  for  further 
studies  to  show  whether  estradiol  insensitivity  occurs  in  clinical 
ER*  LN  metastases  and  whether  other  metastatic  microenviron¬ 
ments  differentially  modulate  estradiol  sensitivity  and  to  define 
estradiol-regulated  genes  uniquely  influenced  by  the  microenvi¬ 
ronment  of  tumor  cells.  This  approach  may  identify  key  genes  that 
can  be  targeted  therapeutically  at  specific  metastatic  sites. 

In  summary,  we  describe  a  xenograft  model  of  EWD  after  LN 
metastasis  and  show  that  estradiol  regulation  of  genes  is  different 
in  ER*  primary  MGTs  compared  with  their  matched  ER*  LN 
metastases.  Through  global  gene  expression  profiling  and  immu- 
nohistochemical  analyses,  we  identify  genes  that  maintain,  lose,  or 
are  regulated  in  the  opposite  direction  by  estradiol  in  tumors  and 
LN  metastases.  These  results  suggest  that  current  ER-targeted 
therapies  may  differentially  affect  tumor  cells  depending  on  the 
microenvironment  in  which  the  cells  reside.  If  so,  our  studies 
provide  possible  new  therapeutic  targets,  which,  in  combination 
with  standard  hormonal  therapies,  might  improve  responses  of 
advanced  breast  cancers. 
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Abstract  Genome-wide  expression  profiling  has  expe¬ 
dited  our  molecular  understanding  of  the  different  suhtypes 
of  breast  cancers,  as  well  as  defined  fhe  differences  among 
genes  expressed  in  primary  tumors  and  their  metastases. 
Laser-capture  microdissection  (LCM)  coupled  to  gene 
expression  analysis  allows  us  to  understand  how  specific 
cell  types  contribute  to  the  total  cancer  gene  expression 
signature.  Expression  profiling  was  used  to  define  genes  that 
contribute  to  breast  cancer  spread  into  and/or  growth  within 
draining  lymph  nodes  (LN).  Whole  tumor  xenografts  and 
their  matched  whole  LN  metastases  were  compared  to  LCM 
captured  cancer  cells  from  the  same  tumors  and  matched 
LN  metastases.  One-thousand  nine-hundred  thirty  genes 
were  identified  by  the  whole  organ  method  alone,  and  1,281 
genes  by  the  LCM  method  alone.  However,  less  than  1% 
(30  genes)  of  genes  that  changed  between  tumors  and  LN 
metastases  were  common  to  both  methods.  Several  of  these 
genes  have  previously  been  implicated  in  cancer  aggres¬ 
siveness.  Our  data  show  that  whole-organ  and  LCM  based 
gene  expression  profiling  yield  distinctly  different  lists  of 
metastasis-promoting  genes.  Contamination  of  the  tumor 
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cells,  and  cross  reactivity  of  mouse  RNA  to  human-specific 
chips  may  explain  Ihese  differences,  and  suggests  that 
LCM-derived  data  may  be  more  accurate. 

Keywords  Breast  cancer  •  Estrogen  receptors  • 

Lymph  node  •  Lymphatic  vessel  •  Microarray  •  Metastasis  • 
Progesterone  receptors  •  Tumor  xenografts  •  ZsGreen 

Abbreviations 

CK18  Cytokeratin  18 

ER  Estrogen  receptor 

H  &  E  Hematoxylin  and  eosin 

LCM  Laser-capture  microdissection 
LN  Lymph  node 

Introduction 

Breast  cancer  metastasis  kills  over  400,000  worldwide  each 
year  [1].  At  diagnosis,  30-50%  of  primary  breast  cancers 
have  disseminated  to  other  locations  including  the  draining 
lymph  nodes  (LN).  Assessments  of  tumor  aggressiveness 
as  predictors  of  overall  survival  were  historically  confined 
fo  histological  evaluations  by  trained  pathologists.  More 
recently,  molecular  approaches  have  been  shown  to  be 
superior  to  pathological  staging  to  determine  patient 
prognosis  and  treatment  options  [1-5]. 

Through  global  gene  expression  profiling,  breast  cancers 
have  been  shown  to  be  derived  from  at  least  two  func¬ 
tionally  distinct  cell  types,  basal  cells  and  luminal 
epithelial  cells;  each  responds  differently  to  therapeutics 
[6].  From  these  genome- wide  analyses  new  hypotheses 
have  been  proposed,  which  suggest  that  the  ability  of  a 
breast  cancer  cell  to  metastasize  is  intrinsically  encoded 
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and  linked  to  the  same  genetic  alterations  that  initiated  the 
growth  of  the  primary  tumor  [7,  8].  Alternatively,  other 
data  indicate  that  only  certain  ‘rare  cells’  within  a  breast 
tumor  harbor  the  potential  for  tumor  growth  and  metastasis 
[9].  It  is  possible  that  both  views  of  progression  are  correct 
and  that  tumors  that  generate  metastatic  cells  differ  only 
slightly  in  their  transcriptome  from  the  bulk  of  the  cells 
that  comprise  the  primary  tumor. 

A  large  number  of  studies  have  been  conducted  to 
identify  the  genes  that  are  expressed  differently  in  primary 
tumors  compared  to  their  metastases,  and  that  may  there¬ 
fore  contribute  to  increased  aggressiveness.  Interestingly, 
the  published  ‘metastatic  gene  prohles’  are  not  composed 
of  the  same  gene  lists  [10].  Among  other  things,  these 
differences  among  studies  have  been  attributed  to  tissue 
sampling  variations  or  to  different  statistical  analyses  [10]. 

We  have  been  engaged  in  identifying  genes  that  may  be 
responsible  for  the  spread  and/or  growth  of  estrogen 
receptor  (ER)  positive  breast  cancers  to  the  LNs.  To  this 
end  we  determined  if  the  two  commonly  used  methods  to 
isolate  RNA — whole-organ  isolation  and  laser-capture 
microdissection  (LCM) — identify  similar  sets  of  genes  in 
an  experimental  model.  We  hnd  that  using  the  same  tumor 
lines  and  statistical  analyses,  gene  expression  profiles  from 
whole-organs  are  completely  different  from  the  profiles 
using  LCM  isolated  cells,  with  only  minor  overlap  between 
them.  Some  of  the  overlapping  genes  have  been  previously 
shown  to  be  involved  in  cancer  aggressiveness.  We  con¬ 
clude  that  the  methods  used  to  isolate  RNA  can  markedly 
skew  the  results  obtained. 


Materials  and  methods 

Fluorescent  xenograft  model  of  ER  +  LN  metastasis 

ERh-  MCF-7  human  breast  cancer  cells  were  originally 
purchased  from  the  American  Type  Culture  Collection 
(Manassas,  VA).  The  generation  of  fluorescent  MCF-7 
cells  has  been  previously  described  [11].  In  brief,  ZsGreen 
retroviral  particles  (Clontech,  Mountain  View,  CA)  were 
isolated  from  PT-67  packaging  cells  (Clontech),  Altered, 
and  overlaid  onto  MCF-7  cells.  Cells  were  serially  trans¬ 
duced  and  then  subjected  to  G418  selection,  followed  by 
aseptic  fluorescence  associated  cell  sorting  to  isolate  a 
homogeneously  bright  green  subpopulation. 

All  animal  procedures  were  performed  under  a  protocol 
approved  by  the  University  of  Colorado  Institutional  Ani¬ 
mal  Care  and  Use  Committee.  Ovariectomized  athymic  nu/ 
nu  mice  were  obtained  from  Harlan  Sprague-Dawley 
(Indianapolis,  IN)  at  5-6  weeks  of  age.  To  establish 
tumors,  mice  were  anesthetized  with  Avertin  and  injected 
into  the  opening  of  the  lactiferous  duct  of  the  abdominal 


mammary  glands  with  1  million  ZsGreen-expressing 
MCF-7  cells  in  100  microliters  of  100%  Matrigel  (Becton 
Dickinson,  Franklin  Lakes,  NJ).  Mice  were  also  implanted 
with  silastic  pellets  containing  17/1-estradiol  (2  mg  +  8  mg 
cellulose)  as  previously  described  [12].  Fluorescent  whole- 
body  imaging  (Illumatool  9900,  Lightools  Research,  En¬ 
cinitas,  CA)  was  done  weekly  to  determine  when  axillary 
LN  metastases  arose.  Mice  were  euthanized  by  CO2 
asphyxiation  and  fluorescent  intravital  optical  imaging  was 
done  by  coupling  the  Illumatool  with  an  Olympus  (Mel¬ 
ville,  NY)  SZ-61  dissecting  microscope  and  Olympus  C- 
5050  digital  camera.  Fluorescent  tumors  in  mammary 
glands  and  matched  LN  metastases  were  removed. 

Whole-organ  RNA  isolation  and  expression  profiling 

Tumors  in  mammary  glands  and  LNs  containing  metasta¬ 
ses  were  extracted  and  frozen  at  -80°C.  RNA  was  prepared 
with  Trizol  according  to  the  manufacturer’s  specifications 
(Invitrogen,  Carlsbad,  CA).  Samples  were  in  vitro  tran¬ 
scribed,  biotinylated,  hybridized  to  Affymetrix  HG- 
U133H-2  chips  that  interrogate  47,000  transcripts  with 
54,675  probesets  (Santa  Clara,  CA),  and  scanned  on  an 
Agilent  bioanalyzer  2100  (Santa  Clara,  CA). 

Laser-capture  microdissection,  RNA  isolation  and 
expression  profiling 

Mammary  gland  tumors  and  LNs  containing  metastases 
were  harvested,  placed  in  a  Tissue  Tek  cryomold  (FMS, 
Hatfield,  PA),  covered  with  Neg50  frozen  section  medium 
(Richard-Allan  Scientific,  Kalamazoo,  MI),  and  frozen  in 
isopentane  (Sigma,  St.  Louis,  MO)  cooled  by  liquid 
nitrogen.  Eight  micron  sections  were  cut  through  the  entire 
tumor  or  LN  with  a  cryotome  set  at  -20°C.  Each  section 
was  collected  on  an  uncharged  slide  (Fisher),  and  every 
10th  slide  was  stained  for  hematoxylin  and  eosin  (H  &  F) 
or  processed  for  cytokeratin  18  (CK18)  immunofluores¬ 
cence  to  confirm  location  of  the  cancer  cells.  Serial 
sections  of  slides  containing  cancer  cells  were  processed 
through  75%  ethanol,  water,  75,  95,  and  100%  ethanol  for 
30  s  each,  and  finally  dehydrated  with  xylene  for  1  min 
before  being  immediately  subjected  to  laser-capture 
microdissection  using  an  Arcturus  Autopix  (Sunnyvale, 
CA).  Two-thousand  cells  for  each  sample  were  extracted 
within  30  min  of  fixation  and  frozen  at  -70°C  until  further 
processing.  Suggested  protocols  from  Arcturus  were  fol¬ 
lowed  for  RNA  extraction  and  two  rounds  of  amplification 
using  the  Picopure  RNA  isolation  kit  and  RiboAmp  HS  kit 
(Arcturus),  respectively.  Samples  were  in  vitro  transcribed, 
biotinylated,  hybridized  to  Affymetrix  (Santa  Clara,  CA) 
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human  X3P  chips  that  interrogate  47,000  transcripts  with 
61,358  probesets,  and  scanned  on  an  Agilent  (Santa  Clara, 
CA)  bioanalyzer  2100. 

Statistical  analyses 

The  data  from  four  whole-organ  chips  (two  pairs  of  tumors 
and  their  matched  LNs),  and  four  LCM  chips  (two  pairs  of 
tumors  and  their  matched  LNs)  were  imported  into  Partek 
(Genomics  Suite  software  St.  Louis,  MO)  and  RMA  nor¬ 
malized.  For  comparative  analyses  between  the  whole- 
organ  and  LCM  data,  matching  probesets  from  the  two 
different  chip  types  were  identihed.  Of  the  54,675  U133-I-2 
probesets,  3,487  were  dropped  from  the  analysis  because 
there  was  no  X3P  match,  but  7,301  probesets  were  used 
multiple  times.  The  result  was  that  all  but  62  of  the  61,358 
X3P  probesets  were  matched  to  a  U133-I-2  probeset.  The 
expression  values  for  each  organ  corresponding  to  either 
the  U133H-2  or  X3P  probeset  were  exported  into  Microsoft 
Excel  (Redmond,  WA).  An  array-comparison  datasheet 
that  matched  the  U 133-1-2  probesets  to  the  X3P  probesets 
was  obtained  from  Affymetrix.  Using  this  program  the 
RMA  normalized  expression  values  from  Partek  for  the 
four  U133-I-2  chips  were  matched  to  the  corresponding  X3P 
probesets  of  the  four  X3P  chips.  This  resulted  in  normal¬ 
ized  expression  levels  from  all  eight  chips.  The  chip  data 
were  also  imported  into  Genespring  (Agilent)  to  obtain 
present/marginal/absent  calls  for  each  chip.  Probesets  with 
absent  calls  on  all  eight  chips  of  the  same  type  were 
excluded  from  analysis  leaving  36,682  genes  that  were 
present  on  one  of  the  eight  chips  for  analysis.  These  genes 
were  imported  into  Partek  according  to  the  X3P  probeset 
identifiers.  Separate  paired  t-tests  were  performed  on  the 
whole-organ  and  LCM  matched  primary  tumor  and  LN 
pairs.  P-values  <0.05  were  considered  signihcant. 

Gene  networking  analyses 

To  dehne  relationships  among  genes  that  were  differen¬ 
tially  regulated  in  tumors  versus  LN  metastases,  pathway 
analyses  were  performed  using  Ingenuity  software  (Red¬ 
wood  City,  CA).  Probeset  identification  numbers  from 
transcripts  with  increased  or  decreased  expression  in  LN 
metastases  were  used  for  gene  identihcation. 

Results 

Tissue  acquisition 

Estrogen-dependent  MCF-7  tumor  growth  is  required  for 
LN  metastasis  [11].  For  these  studies  MCF7  +  ZsGreen 


cells  were  grown  in  the  mammary  gland  of  ovariectomized 
immunocompromised  mice  that  had  been  supplemented 
with  an  estradiol  releasing  pellet.  Identification  of  axillary 
LN  metastases  was  achieved  with  fluorescent  whole  body 
imaging  (Fig.  la).  At  12  weeks  intact  mammary  gland 
tumors  (Fig.  lb)  and  their  matched  LN  metastases  (Fig.  Ic) 
were  harvested  and  frozen  as  intact  whole-organs  or  cut 
into  8-micron  frozen  sections.  LNs  and  tumors  were  then 
probed  with  a  breast  cancer  cell  marker,  cytokeratin  18,  to 
map  cancer  cell  locations  (Fig.  Id).  Adjacent  sections  were 
then  used  for  LCM  (Fig.  le-g).  The  LN  metastases  that 
were  used  for  the  intact  whole-organ  analysis  contained 
approximately  equal  amounts  of  cancer  cells  as  assessed  by 
ZsGreen  fluorescence,  CK18  immunohistochemistry,  and 
H&E  staining. 

Gene  expression  profiling 

To  compare  gene  expression  signatures  in  whole  tumors 
and  matched  intact  LN  metastases,  a  paired  t-test  was 
performed.  With  this  approach  1 ,930  genes  were  identihed 
as  being  differentially  expressed  in  tumors  compared  to  LN 
metastases,  with  1,213  decreased  and  717  increased  in  the 
metastases  compared  to  the  tumors  from  which  they  orig¬ 
inated  (Fig.  2a).  For  the  LCM  samples,  the  same 
procedures  were  followed,  which  identihed  368  decreased 
transcripts  and  913  increased  transcripts  in  the  metastases 
compared  to  the  tumors  from  which  they  originated 
(Fig.  2b).  Clearly  the  pattern  of  upregulation  and  down- 
regulation  differs  in  the  two  types  of  analyses. 
Additionally,  the  volcano  plots  show  much  larger  differ¬ 
ences  in  ‘fold-changes’  between  tumors  and  LN  metastases 
that  were  laser-captured  compared  to  the  whole-organ 
samples. 

Genes  that  change  in  whole-organ  compared  to  LCM 

When  whole-organ  and  LCM  paired  t-test  results  were 
combined,  3,181  transcripts  were  found  to  have  signih- 
cantly  increased  or  decreased  expression  levels  in  tumors 
versus  LN  metastases.  Of  the  3,181  however,  only  68  (2%) 
were  identihed  as  changing  by  both  analyses  (Fig.  3a),  and 
when  the  direction  of  change  was  calculated  for  these  68 
transcripts;  38  were  found  to  change  in  opposite  directions 
by  the  two  analyses  (Fig.  3a).  The  remaining  30  transcripts 
(~  1%  of  all  genes)  that  changed  in  the  same  direction  in 
both  the  whole-organ  and  LCM  analyses  are  listed  in 
Table  1. 

Results  from  the  paired  t-tests  were  compared  to  identify 
all  transcripts  that  were  differentially  regulated  in  either 
whole-organ  or  LCM  analyses.  Using  both  data  sets,  1,570 
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Fig.  1  Isolation  of  tumors  and 
LN  metastases  used  for  “whole- 
organ”  or  laser-capture 
microdissection  followed  by 
gene  expression  profiling,  (a) 
Fluorescent  whole  body  image 
of  an  immunocompromised 
mouse  with  bilateral 
MCF7  +  ZsGreen  tumors  and 
an  axillary  LN  metastasis,  (b) 
ex  vivo  fluorescent  image  of  a 
MCF7  4-  ZsGreen  tumor,  (c) 

LN  metastasis,  (d)  cytokeratin 
18  immunohistochemistry  of  an 
axillary  LN  metastasis,  CK18 
red,  DAPI,  blue,  (e) 
hematoxylin  and  eosin  stain  of  a 
LN  identifying  the  metastatic 
MCF7  4-  ZsGreen  cells  (*),  (f) 
combined  images  of  a  serial 
section  of  d  with  the  cancer 
cells  marked  for  LCM,  (g)  Post- 
LCM  image  of  captured  MCF7 
cells  from  e 


transcripts  were  found  to  be  downregulated  in  LN  metas¬ 
tases  (Fig.  3b).  Of  these,  1,202  were  unique  to  the  whole- 
organ  analysis,  357  were  unique  to  LCM  analysis,  and  1 1 
genes  exhibited  decreased  expression  with  either  method. 
A  similar  trend  was  observed  among  1,61 1  genes  that  were 
upregulated  in  LN  metastases,  with  698  being  whole-organ 
specific,  894  being  LCM  specific,  and  19  common  to  the 
two  analyses  (Fig.  3c). 

Ingenuity  pathway  analyses  were  performed  to  deter¬ 
mine  if  the  large  differences  between  the  two  types  of 
analyses  could  nevertheless  be  due  to  changes  among 
similar  functional  networks.  The  whole-organ  data  identi¬ 
fied  pathways  encoding  GM-CSF  signaling,  insulin 
receptor  signaling,  amyloid  processing,  Jak/Stat  signaling, 
and  sonic  hedgehog  signaling  that  were  the  most  different 
in  the  LN  metastases.  The  LCM  data  identified  completely 


different  pathways  that  had  changed  the  most  in  LN 
metastases  including  PI3K/AKT  signaling,  PTEN  signal¬ 
ing,  synthesis  and  degradation  of  ketone  bodies,  oxidative 
phosphorylation,  and  keratin  sulfate  biosynthesis.  Clearly 
both  the  individual  genes  and  the  broader  pathways  that  are 
altered  in  tumors  versus  LN  metastases  are  different 
depending  on  the  assays  that  are  performed. 

Discussion 

Numerous  studies  from  independent  labs  have  shown  that 
there  are  gene  expression  changes  in  cancer  cells  within 
primary  tumors  compared  to  cells  that  seed  to  different 
metastatic  microenvironments  [8,  13].  These  studies  sug¬ 
gest  that  the  propensity  for  development  of  metastasis  is 
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a  Change  in  LN  metastases  gene  expression 
(whole-tissue) 

Decreased  Increased 


-4  -2  N/C  2  4 

Fold  change  LN  metastases/tumors 


b  Change  in  LN  metastases  gene  expression 

(laser-captured) 

Decreased  Increased 


-4  -2  N/C  2  4 

Fold  change  LN  metastases/tumors 


Fig.  2  Identification  of  the  magnitude  of  gene  expression  changes  in 
tumors  and  LN  metastases.  (a)  Volcano  plot  showing  the  1,930  genes 
that  were  significantly  different  between  tumors  and  LN  metastases  in 
the  whole-organ  analysis,  (b)  volcano  plot  showing  the  1,281  genes 
that  were  significantly  different  between  tumors  and  LN  metastases 
by  the  LCM  analysis.  For  both  plots  the  magnitude  of  fold  change  is 
shown  on  the  x-axis  and  the  significance  is  shown  on  the  y-axis 

due  to  genes  that  are  altered  during  the  initial  tumorigenic 
process.  The  genetic  changes  in  metastases  have  been 
attributed  to  differences  required  for  successful  coloniza¬ 
tion  and  growth  within  vastly  different  metastatic 


a  Genes  significantly  different  in  LN 
metastases  by  both  analyses 


Decreased  Opposite  Increased 


b  Genes  with  decreased 

expression  in  LN  metastases 


Whole  tissue  LCM 


C  Genes  with  increased 

expression  in  LN  metastases 


Whole  tissue  LCM 


Fig.  3  Determining  the  overlap  of  the  genes  that  changed  in  the 
whole-organ  and  LCM  analyses.  Significantly  changing  genes  have  a 
decrease  or  increase  in  the  direction  of  the  fold  change  without 
specifying  a  cut-off  value,  (a)  Categorization  of  the  68  genes  that 
were  changed  in  both  whole-organ  and  LCM  analyses.  Shown  are 
genes  that  decreased  or  increased  in  LN  metastases  with  both  analyses 
or  those  whose  expression  changed  in  the  opposite  direction  in  each 
analysis,  (b)  genes  with  decreased  expression  in  LN  metastases,  (c) 
genes  with  increased  expression  in  LN  metastases 


microenvironments  [14].  Whether  the  cells  that  are  capable 
of  metastasis  are  ‘rare  cells’  or  ‘cancer  stem-cells’,  or 
reflect  any  cell  within  the  tumor,  remains  under  debate  [8, 
9].  Regardless,  determining  the  genes  that  are  differentially 
expressed  between  cancer  cells  growing  within  the  breast 
compared  to  those  growing  in  other  sites  is  crucial  to 
understanding  the  mechanisms  that  allow  cells  to  seed  and 
grow  at  distant  sites.  Identifying  these  genes  and  rendering 
them  therapeutic  targets,  is  an  essential  first  step  towards 
inhibiting  lethal  metastases. 
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Table  1  Genes  that  increased  or 

decreased  in  LN  metastases  by  both  analyses 

X3P  PSID 

UniGene  ID 

Gene  symbol 

Whole-tissue  FC(LN/MGT) 

LCM  EC(LN/MGT) 

Genes  with  increased  expression 

in  LN  metastases 

in  both  analyses 

219529_3p_at 

Hs.64746 

CLIC3 

1.1114 

1.0972 

Hs2.429028.1.Sl_3p^at 

Hs.632184 

CARHSPl 

1.0548 

1.0498 

gll545730_3p_at 

Hs.  1 12569 

GAN 

1.0492 

1.1879 

g4557696_3p_a_at 

Hs.99936 

KRTIO 

1.0431 

1.2789 

208549_3p_x_at 

Hs.591909 

PTMA 

1.0376 

1.1448 

g4507574_3p_a_at 

Hs.279594 

TNFRSFIA 

1.0373 

1.1685 

g388890_3p_a_at 

Hs.  89643 

TKT 

1.0285 

1.2193 

g9295 1 9  l_3p_s_at 

Hs.584881 

C3orf28 

1.0242 

1.3023 

gl  2804874 Jp_a_^at 

Hs.226390 

RRM2 

1.0222 

1.3370 

g4505038_3p_at 

Hs.1116 

LTBR 

1.0181 

1.1803 

Hs.l21017.0.Sl_3p_at 

Hs.121017 

HIST1H2AE 

1.0171 

1.3161 

g7959918_3p_a_at 

Hs.381126 

RPS14 

1.0153 

1.2751 

Hs.27931.0.Sl_3p_at 

Hs.592295 

GNPNATl 

1.0083 

1.1926 

g7705770_^3p_at 

Hs.  19385 

ABHD5 

1.0062 

1.3309 

Hs.27239.0.Sl_3p_at 

Hs.27239 

ZDHHC5 

1.0057 

1.1618 

201217_3p_x_at 

Hs.  1 19598 

RPL3 

1.0050 

1.1517 

211073Jp_x_at 

Hs.  1 19598 

RPL3 

1.0043 

1.1548 

gl3623708_3p_s_at 

Hs.  1 19598 

RPL3 

1.0042 

1.1483 

g450 1 976_3p_a_at 

Hs.512815 

AP3D1 

1.0028 

1.2395 

Genes  with  decreased  expression 

in  LN  metastases 

in  both  analyses 

Hs.74034.1.SC3p_at 

Hs.74034 

CAVl 

-1.14927 

-1.21563 

Hs.24812.0.Sl_3p_at 

Hs.472027 

CDS2 

-1.01957 

-1.06974 

Hs.l72329.1.Sl_3p_at 

Hs.485628 

ZNF451 

-1.01761 

-1.06306 

Hs.l61962.0.Sl_3p_at 

Hs.370510 

IGSF4 

-1.03359 

-1.01245 

Hs.l45612.0.S2„3p_at 

Hs.  190622 

DDX58 

-1.05818 

-1.26471 

Hs.l06440.0.Al_3p_at 

- 

- 

-1.07679 

-1.05797 

g9951924_3p_at 

Hs.210995 

CA12 

-1.08361 

-1.05245 

g7661679_3p_a_at 

Hs.455336 

NELE 

-1.01845 

-1.04308 

g4929370_3p_x_at 

Hs.465885 

ILF3 

-1.0444 

-1.0716 

gl2052731 _3p_x_at 

Hs.351906 

ZNF506 

-1.04132 

-1.0935 

gl2052731_3p_at 

Hs.351906 

ZNF506 

-1.04186 

-1.07409 

Shown  are  the  probeset  identification  code  (PSID),  the  UniGene  identification  code,  the  gene  symbol  and  the  RMA-normalized  fold  changes 
from  the  whole-organ  and  LCM  analyses 


For  the  whole-organ  analysis,  RNA  was  isolated  from 
entire  tumors,  or  from  whole  LNs  containing  similar 
amounts  of  cancer  cells  within  them.  For  LCM,  RNA  was 
isolated  from  2,000  pure  cancer  cells  of  tumors  and  LN 
metastases  (Fig.  1).  To  compare  these  sets  we  chose  an 
approach  that  is  commonly  used  to  study  large  numbers  of 
gene  expression  changes;  namely,  global  gene  expression 
profiling  with  Affymetrix  gene  chips.  The  chips  used  are 
designed  to  theoretically  hybridize  only  with  human  RNA, 
and  therefore  both  the  whole-organ  and  LCM  approaches 
should  inform  us  only  of  the  gene  expression  changes  that 
occur  in  the  human  cancer  cells  within  the  mouse  host 
sites.  It  was  therefore  surprising  to  find  greater  ‘fold 


changes’  in  LN-specific  genes  assessed  by  LCM  (Fig.  2a 
vs.  Fig.  2b). 

Since  the  gene  chips  used  for  the  two  methods  interro¬ 
gated  the  same  transcripts,  with  all  but  62  of  the  61,358 
X3P  probesets  having  a  matching  U133+2  chip  probeset,  it 
was  surprising  that  very  few  genes  changed  similarly 
between  the  two  analyses  (Fig.  3).  Furthermore,  of  the  68 
transcripts  that  were  significantly  different  in  LNs  by  both 
analyses,  38  were  regulated  in  the  opposite  direction,  being 
increased  in  the  metastases  by  one  analysis,  and  decreased 
by  the  other  (Fig.  3a).  These  results  suggest  that  the  Af¬ 
fymetrix  chips  are  binding  not  only  human  RNA  but  also 
mouse  RNA.  This  aberrant  hybridization  may  be 
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responsible  not  only  for  the  gene  differences  reported  by 
the  two  methods  but  also  for  the  lack  of  overlap  between 
the  functional  pathways  identified. 

No  other  reports  have  determined  how  sample  collection 
methods  influence  LN  metastasis  gene  signatures.  How¬ 
ever,  one  study  [15]  compared  28  breast  tumors,  17  ER 
positive  versus  11  ER  negative,  and  identifled  genes  dif¬ 
ferently  expressed  between  these  two  cancer  types  using 
both  LCM  and  bulk-tumor  gene  expression  profiling. 
Interestingly,  of  all  the  genes  found  to  be  significantly 
different  between  ER  positive  and  ER  negative  tumors,  less 
than  one-third  were  common  to  both  analyses.  These 
findings  are  interesting  because  the  same  human  tumors 
were  used  for  comparison,  yet  the  ER  signature  changes 
were  dependent  on  the  type  of  analysis  conducted. 

We  suggest  two  possible  explanations  for  the  differ¬ 
ences  between  our  whole  organ  and  LCM-generated  data. 
(1)  It  is  possible  that  multiple  unique  transcripts  are  altered 
by  the  amplification  procedure  required  for  LCM.  How¬ 
ever,  a  previous  study  showed  that  amplification  from  the 
3'  end  of  transcripts  does  not  have  a  major  impact  on 
overall  gene  profiles  with  Affymetrix  GeneChips  because 
the  probe  sets  on  the  arrays  are  designed  using  3'  end 
sequences  [16].  Additionally,  another  report  found  no 
significant  differences  in  the  expression  intensity  of  21 
housekeeping  genes  from  matched  LCM  and  bulk-tumor 
gene  expression  data  [15].  (2)  An  alternative  hypothesis 
posits  that  Affymetrix  chips  can  bind  not  only  human 
RNA,  but  also  mouse  RNA,  thus  altering  the  chip  readout 
from  whole  organ  analyses.  Naef  et  al.  [17]  recently  tested 
this  idea  and  found  that  indeed,  mouse  RNA  hybridized  to 
the  “human-specific”  gene  chips.  Interestingly,  we  found 
that  on  average  16.4%  more  probesets  were  present  on 
whole-organ  samples  compared  to  LCM  samples,  sug¬ 
gesting  that  in  xenografts,  mouse  RNA  was  contaminating 
the  whole-organ  analysis.  The  studies  by  Naef  et  al.  sug¬ 
gest  not  only  the  possibility  of  cross-species  contamination, 
but  the  possible  skewing  of  data  due  to  intraspecies  con¬ 
tamination  of  tumor  cells  by  normal  epithelial,  stromal, 
vascular,  and  immune  tissue  elements  within  an  organ. 
Thus  it  may  be  critical  to  use  LCM  to  determine  how  each 
tissue  component  of  an  organ  contributes  to  the  overall 
gene  signature.  Based  on  this,  we  suspect  that  our  LCM 
data  represent  a  more  accurate  gene  expression  comparison 
between  tumors  and  their  matched  LN  metastases. 

These  studies  nevertheless  Identified  30  transcripts  that 
were  increased  or  decreased  in  LN  metastases  by  both 
analyses  (Table  1).  These  genes  may  represent  the  best 
targets  for  further  studies  to  identify  the  mechanisms  that 
contribute  to  the  spread  or  growth  of  breast  cancer  cells 
into  and  within  LNs.  Prothymosin  alpha  (PTMA)  was 
increased  in  the  LN  metastases  and  has  previously  been 
shown  to  increase  proliferation  rates  of  breast  cancer  cells 


[18].  This  finding  is  interesting  in  view  of  our  previous 
findings  of  increased  proliferation  rates  in  LN  metastases 
[11].  Tumor  necrosis  factor  receptor  superfamily  member 
lA  (TNERSEIA)  was  also  increased  in  the  LN  metasta¬ 
ses.  Its  expression  has  previously  been  shown  to  decrease 
apoptosis  in  MCF-7  cells  treated  with  human  TNE  alpha 
protein  [19].  Ribonucleotide  reductase  M2  polypeptide 
(RRM2)  and  lymphotoxin  beta  receptor  (LTBR)  were 
both  increased  in  the  LN  metastases,  and  both  increase 
cancer  cell  survival  [20,  21].  Interestingly,  RRM2  is  tar¬ 
geted  by  gemcitabine  and  hydroxyurea  in  lung  and  skin 
cancers,  respectively  [22].  A  transcript  whose  expression 
was  decreased  in  LN  metastases  and  has  been  studied  for 
its  effects  on  cancer  progression  is  caveolin-1  (CAVl). 
This  scaffolding  protein  is  the  main  component  of  the 
caveolae  in  plasma  membranes  of  most  cell  types,  and  its 
expression  inhibits  proliferation  in  human  breast  cancer 
cells  [23].  Its  decreased  expression  may  therefore  play  a 
role  in  the  increased  cell  proliferation  seen  in  LN 
metastases  [11,  24]. 

In  summary,  this  report  evaluated  whole-organ  versus 
LCM  for  gene  expression  profiling,  in  order  to  identify 
therapeutic  targets  that  may  reduce  breast  cancer  spread  to 
and  growth  in  LNs.  However,  when  comparing  these  two 
experimental  approaches,  we  found  almost  no  overlap 
among  the  transcripts  and  functional  pathways  identified. 
The  30  transcripts  that  were  common  to  both  analyses  may 
be  the  most  promising  for  future  therapeutic  targeting. 
However,  our  data  clearly  demonstrate  that  the  methods 
used  to  prepare  tumors  for  gene  expression  profiling  can 
markedly  skew  the  results  obtained. 
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